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,  20,  Abstract  (continued) 

are  described  and  compared  to  data  from  the  recovered  experimental  model.  The 
tunnel  closure  predictions  and  the  predicted  locations  of  block  motion  are 
in  reasonable  agreement  with  some  experimental  data. 

Other  aspects  of  the  study  included: 

•  development  and  application  of  an  analytical  technique  for 
predicting  joint  activation  conditions  (i.e.,  block  motions) 
for  specified  joint  orientations  and  locations  near  the 
tunnel  surface, 

•  calculations  and  comparisons  of  tunnel  response  to  dynamic 
and  quasi-static  stress  loads,  and 

•  investigation  of  the  effects  of  joint  friction  on  tunnel 
damage. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 


1 . 1  BACKGROUND 

All  rock  masses  contain  large  numbers  of  planar  discontinu¬ 
ities  or  joints.  These  joints  often  react  to  an  applied  stress 
quite  differently  than  the  surrounding  material  and  there  is  both 
experimental  (Ref.  1)  and  theoretical  (Refs.  2,  3)  evidence  to 
suggest  that  sliding  motion  across  these  joints  should  be  one  of 
the  important  considerations  involved  in  designing  underground 
structures  which  are  to  survive  nuclear  and/or  conventional 
explosions. 

In  the  recent  DIABLO  HAWK  underground  nuclear  test,  SRI 
International  fielded  several  laboratory-scale  structures  in 
jointed  rock  simulant  (Ref.  4).  Figure  1  shows  the  general  model 
designs.  The  jointed  rock  models  are  18"  and  30"  in  diameter, 
with  3"  an^  5"  diameter  tunnels,  respectively.  The  SRI  small- 
scale  jointed  structures  experiment  provides  an  opportunity  for 
examining  the  influence  of  highly  idealized  jointed  arrays  on 
tunnel  closure  and  damage  characteristics. 

There  are  both  continuum  and  explicit  approaches  for  pre¬ 
dicting  tunnel  closure  and  damage  in  jointed  media.  In  the 
continuum  approach,  the  overall  mechanical  properties  of  the 
field  are  synthesized  from  the  properties  of  the  intact  rock  and 
the  properties  and  arrangement  of  the  joints;  for  example,  stiff¬ 
nesses  and  strengths  used  in  the  continuum  model  are  lower  than 
the  intact  rock  properties.  In  the  explicit  approach,  the  prop¬ 
erties  and  arrangement  of  individual  blocks  and  joints  are 
individually  and  explicitly  modeled  and  used  to  calculate  tunnel 
response  to  loading.  Both  of  these  treatments  have  important 
limitations:  The  rationale  for  choosing  a  set  of  degraded 


s> 
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propLi ties  in  the  continuum  approach  is  usually  hard  to  justify, 
and  this  approach  precludes  discontinuities  in  displacements 
across  joints,  which  may  be  a  critical  omission.  The  properties 
of  joints  in  the  explicit  approach  are  hard  to  measure,  and  the 
arrangement  of  joints  in  a  real  field  around  a  real  tunnel  can 
only  be  estimated. 


1..’.  OBJECTIVE-:  AND  APPROACH 


The  >b  native  of  this  study  is  to  examine  the  influence  of 
the  joints  on  tunnel  closure  and  damage  using  numerical  techniques 
which  explicitly  treat  the  idealized  joints  in  the  SRI  structure. 

The  approach  involves  three  related  aspects: 

1.  Joint  Activation. 

An  analytical  closed  form  solution  is  derived  for  the 
initiation  or  activation  of  relative  motion  across  a 
joint.  The  joint  activation  depends  on  (a)  the  rela¬ 
tive  orientation  of  the  joint  and  the  tunnel  surface 
intersected  by  the  joint,  (b)  the  magnitude  and  propa- 
uation  direction  of  the  stress  wave  encompassing  the 
tunnel,  and  (e)  the  shear  stress  on  the  joint  which,  must 
exceed  a  critical  value  which  depends  on  tne  normal 
stress  across  the  joint  as  well  as  on  surface  roughness. 

In  Section  2 ,  the  joint  activation  solution  is  derived 
anti  applied  to  the  SRI  small-scale  models  to  estimate 
(a)  the  areas  on  tne  tunnel  surface  where  joint 
activation  is  likely  to  occur,  (b)  the  stress  levels 
which  are  necessary  to  activate  joints  of  various 
orientations,  and  (c)  the  sensitivity  of  joint 
activation  to  material  properties. 
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2 .  Dynamic  Tunnel  Response  Cases . 


In  the  SRI  small-scale  models  with  tunnel  diameters  of 
3"  ,  the  transit  time  across  this  characteristic  dimen¬ 
sion  is  only  22  usee  in  the  16A  rock  simulant  material. 
Thus ,  since  the  rise  times  to  peak  stresses  of  0.8  kbar 
in  tiie  DIABLO  HAWK  event  are  5-10  msec,  those  small- 
scale  structures  are  actually  loaded  (and  unloaded)  in 
a  quasi-static  fashion.  Full-scale  tunnels  may  be 
loaded  more  dynamically,  therefore  the  effects  of 
dynamic  loading  on  small-scale  tunnels  buried  in  jointed 
media  are  analyzed  using  six  plane  strain  computer  code 
simulations  with  dynamic  loading  times  of  100-400  usee. 
These  cases  also  vary  the  assumed  joint  and  tensile 
strength  properties  of  the  rock  simulant.  The  WAVE-L 
Lagrangian  code  is  used;  the  basic  formulation  is 
similar  to  the  HEMP  code  (Ref.  3) .  The  numerical 
solutions  and  material  properties  are  described  in 
Section  3. 

3.  Ouas i -_S tat_ic  Tunnel  Response  Prediction  For  a  Small-Scale 
Structure  Loaded  to  0.8  kbar  During  the_ DIABLO  HAWK 
Even t . 

A  numerical  simulation  prediction  of  tunnel  closure  and 
damage  is  made  and  compared  with  a  fielded  small-scale 
SRI  1  ( t  he  ;  u<l  has  a  3"  tunnel  diameter  with  F " 

In*,  quasi  -st  it  i  •  so  i  u  t  i  on  technique  and  the  c  aiculal  lonal- 
r iment  1 1  con;;  ari  sons  are  described  in  Sect  ion  4. 


1.3  SUMMARY  AMD  CONCLUSIONS 
1.3.1  Joint  Activation 

Although  the  response  of  a  tunnel  in  a  jointed  medium  must 

be  solved  numerically,  the  conditions  that  will  initiate  joint 

slipping  can  be  analyzed  in  closed  form  if  we  assume  an  elastic 

media  with  joints,  and  a  far  field  stress  state  of  uniaxial 

strain  with  a  stress  loading  of  o, 

load 

Figure  2  shows  the  loading  stress  required  to  activate 

a  joint  oriented  at  a  =  0,  ‘45°,  and  90°  as  a  function  of  the 
angle  (3)  wiicre  the  joint  intersects  the  tunnel  surface.  Thus 
a  joint  oriented  at  i  =  45°  will  require  a  loading  stress  of 
.06  kbar  (  900  psi)  to  become  activated  if  this  joint  intersects 
tiie  tunnel  surface  at  3  =  100°  (see  the  following  sketch  and 
Figure  2 )  . 


N  Behavior  for  F  +  180° 
\ 

N 
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(  s  i  .1.1 1:  p.ipiM|s  u  r  silo  i  i  i  pii,, 


Note  that  for  each  joint  orientation  (■<),  the  shaded  reqion 

on  Figure  2  shows  the  angles  (.  )  on  the  tunnel  surface  where  joint 

activation  is  possible  as  a  function  of  loading  stress  {  loacj)- 

Note  that  symmetry  implies  equivalent  tunnel  behavior  at  angles 

of  +  180°.  Joints  will  not  slip  outside  of  these  shaded  regions. 

The  tunnel  angles,  .  ,  associated  with  absolute  minimums  in  load- 

inn  stress  (■  ,  . )  nocessarv  for  joint  activation  are  indicated 

■  load 

in  Table  1.  However,  note  that  the  activation  stress  curves 
(Figure  2)  are  relatively  flat  bottomed;  and  therefore,  joint 
activation  can  occur  over  a  fairly  wide  range  in  tunnel  angles 
with  only  a  small  increase  in  loading. 

Table  1.  Tunnel  Angles  (,)  for  Minimum  Loading 

Stress  Lor  Various  Joint  Orientations  ( () . 


Joint 

1  Orientation 

funnel 

Angles  for  Minimum 
oading  Stress* 

j 

l 

ci ci <.  j  roe 

0 

degree 

,  1 

,  .  load 

Kbar 

psi 

07,  113 

.068 

9  86 

:  4  5 

103 

.0  59 

8  5b 

:  ijo 

i 

41,  139 

.102  2 

350 

, . 

> 

.  .  ....  .i 

7 o  ,  103 

.059 

8  5b 

*  Representative  constitutive  properties  are  used 
(fuse  1  ’  1  )  :  Joint  folios  ion  (io)=.04  kbar,  Friction 
Annie  (  : )  ■  11° ,  Material  Poisson's  Ratio  ( v )  =  .  2  5  . 

Figure  >  shows  tlie  areas  on  the  tunnel  surface  where  joint 
ill  iv.it  ion  (for  u  single  joint)  could  occur  for  ,  =  .5  kbar. 
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All  the  joint  orientations  in  the  SRI  experiment  are  shown  for  a 
representative  set  of  joint  and  earth  media  properties.  (Varia¬ 
tions  in  the  joint  and  earth  media  constitutive  properties  are 
illustrated  in  Section  2.) 

1.3.2  Dynamic  Response  Cases 

Six  WAVE-L  numerical  simulations  (Table  3)  were  performed 
involving  100-400  usee  stress  rise  times  to  a  peak  loading  stress 
of  0.8  kbar  on  the  SRI  3"  diameter  models.  The  wave  transit  time 
across  the  3"  tunnel  opening  is  22  |isec,  therefore,  the  dynamic 
tunnel  response  characteristics  could  be  investigated. 

The  tensile  stress  history  near  the  tunnel  crown  is  strongly 
dependent  on  tiie  stress  wave  rise  time  for  times  corresponding 
to  a  I  t'W  transit  times  across  the  tunnel  diameter.  A  stress 
wave  rise  t  i  ::k  •  oi  100  nsec  (-5  transit  times)  caused  a  peak  ten¬ 
sile  at  i  ess  el  almost  1j0  tiers.  In  the  quasi-static  loading  case 
w i t neat  ini  lie  nut,  no  tensile  stress  developed  near  the  crown. 

.  e  a .  ;  ..  i  icain  :  .at  the  material  near  the  tunnel  surface  is 
a  ;s.>  t>  i  it  i  ;  : .  >  ;u.  i  <  lative  displacements*  across  the  joints. 

Wuen  .a  i;  ;  sail  s  across  joints,  the  plates  of  rock  separated 

i  .a  ,  i : : ‘  .  ms  act  as  unsupported  beams  due  to  the  presence 
i  ■ :  t  tie  an  i  i  i.ea  tunnel  .  'a,  result  ina  beam  bend  ing  can  lead  to 

‘  •  • '  a  i  !  •  s:  ;  ,  u.  i  i  r.icture. 

1  .  i.  1  i.tiasi  -  SI  it  ic  Tunnel  Response  Prediction 

■;  lie  pa.  is  i  -  lit  at  ic  solution  technique  is  described  in  Section  4. 
h  s::.e  r  meal  s  i  mu  l  at  ion  pi  e.i  i  ct  ions  using  this  technique  were  made 
;  g  a  !  lelavd  sm.il  1  -sea  1  e  (5"  d  i  a  me  t  o  r  tunnel)  SRI  model.  The 

*  Tiie  leluLivc  i  i  s;  •  1  acemen  t.  or  slipping  across  a  joint,  is  measured 
by  t.  nc  current  distance  between  two  points  which  were  initially 
t. aliening  eaeh  otii.-i  on  opposite  faces  of  the  joint. 


calculation  predictions  were  made  prior  to  excavation  and  examin¬ 
ation  of  the  model. 

Figure  4  shows  a  comparison  of  the  predicted  tunnel  closure 
(as  a  %  of  the  diameter)  versus  tunnel  angle,  .  The  agreement 
is  reasonable  except  near  the  crown-invert  (5  =  0°  and  180°) . 

The  maximum  predicted  closure  is  1.2?.  at  o  =  30°.  The  experi¬ 
mental  data  show  a  maximum  1.6%  closure  at  the  crown-invert 
compared  to  a  0.5%  predicted  value  at  the  crown-invert.  This 
difference  may  be  due  to  the  moderate  sliding  on  the  joints  and 
associated  tensile  failure  and  relaxation  which  occurs  in  the 
calculation,  but  not  in  the  experiment. 

In  the  post-shot  experimental  model,  there  was  deformation 
of  a  thin  (.003  inch)  stainless  steel  tunnel  liner  which  surrounds 
the  foam  filler.  This  liner  may  have  inhibited  sliding  along  the 
joints.  The  observed  liner  deformation  is  associated  with  the 
joints  which  are  predicted  to  slide  the  most;  thus,  some  "block 
motion"  appears  to  have  occurred. 

It  is  concluded  that  tunnel  closure  in  a  well-characterized 
jointed  media  could  be  predicted  using  existing  computer  code 
techniques  involving  explicit  descriptions  of  joint  locations, 
orientations,  and  friction  properties.  Also,  because  of  the 
important  role  joint  (and  fault)  motion  plays  in  determining 
damage  to  tunnel  liners  and  nearby  earth  material,  it  is  felt 
that  continuum  approaches  to  predicting  tunnel  closure  and  damage 
might  not  be  as  useful. 
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SECTION  2 


JOINT  ACTIVATION  ANALYSIS 


The  necessary  conditions  for  a  joint  or  block  to  begin  slip¬ 
ping  depend  on  the  relative  of  the  joint,  the  tunnel 

surface,  and  the  stress  wave,  as  indicated  in  the  following 
sketch : 


\ 

\ 


<  is  the  joint  orienta¬ 
tion  angle 

is  the  tunnel  angle 
(angle  where  joint 
interseets  tunnel) 

is  the  angle  between 
the  joint  and  the 
normal  to  the  tunnel 

s nr iaee ,  .  =  ;  ’ 


Joint  slippage  will  occur  when  the  shear  stress  (i)  on  the  joint 

exceeds  a  critical  value,  i  .  This  value  depends  on  the  normal 

max  1 

stress  (■  )  across  the  joint  as  well  as  on  the  surface  roughness. 
Prior  to  joint  activation,  the  earth  media  around  the  tunnel  is 
assumed  to  respond  elastically  and  statically  to  a  uniaxial  strain 
load  1  ream  a  (guasi-static )  p  Lanai:  wave. 


1  b 


■  i* 


The  key  characteristics  of  interest,  with  respect  to  Joint 

;  '  are  : 


•  the  areas  on  the  tunnel  surface  where  joint 
activation  is  likely  to  occur, 

•  the  stress  levels  which  are  necessary  to 
activate  joints  of  various  orientations, 

•  the  sensitivity  of  joint  activation  to  material 
properties . 


Tne  determination  of  these  characteristics  involves  the 
following  elements: 


1.  The  ti.ii  /*.  field  on  the  unlined 

tunnel  surface  (this  stress  field  is  shown  in 
Figure  5  as  a  ratio  with  the  loading  stress, 
i.e.,  o  /  0  ioacj) 

2.  The  no.J  (i.e.,  we  assume 

i  =  i  +  (tan  4i)u  where  i  and  are 
max  o  n  o  ' 

experimentally  determined  parameters) . 


j  . 


The 


("?oad> 


o .  n . 


■•  •  •  '  ’  -  :c>;  as  a  function  of  joint  orientation 
(  i)  ,  position  on  the  tunnel  ( J )  ,  and  material 
properties . 


The  activation  stress  load  relationship  (derived  in 
Appendix  A)  is 


whe  re 


a 

load 


( 2  i  qcos j)  /  ( sm  ( 2 y-  ;  )  -  sint) 


2  -  11  +  2cos  (2;-.)  ]  4^ 

1 1- v 


(1) 


j  =  | . x  —  .•>  |  (see  sketcn  on  previous  page) 

i  and  ;  are  joint  constitutive  parameters 
o 

is  the  Poisson's  Ratio  for  the  earth  material 
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Note  that  Equation  1  shows  that  r  ,  4> ,  and  v  are  the  con¬ 
stitutive  properties  of  interest  for  joint  activation.  Table  2 
lists  the  numerical  values  of  these  properties  which  were  used 
to  calculate  joint  activation  relationships. 


Table  2 .  Joint  and  Earth  Media  Constitutive 
Properties 


Joint 


Case  Cohesion*,  u 

c 

kbar  psi 


Earth  Media 


Friction  Angle*  Poisson's  Ratio 


1  (degrees) 


*The  joint  cohesion  (i  )  and  friction  angle  (!)  are  related  to 

the  unconfined  compressive  strength  (o  )  of  the  earth  media  by 

2  cos  : 
o 

u  , 

1  -  sin: 

Figure  2  and  Figures  n  to  9  show  the  variations  in  n , 

J  load 

versus  tunnel  angle  (;■■.)  for  the  joint  and  earth  media  property 
variations  indicated  in  Table  2.  Figure  3  and  Figures  10  to  13 
show  the  corresponding  regions  on  the  tunnel  surface  where  joint 
activation  is  possible  for  o  ,  =  .5  kbar. 


joint  Activation  Stress  Loads  tor  Various  Joint  Orientations  (0  as  a  iuiulion  ol 
i  mint*  1  Angle  (■  )  for  Constitutive  Property  Case  1‘2  (  =.US,  ;  -  11"  ,  ,  =  .2j). 


(OEG)  BETA  COEG) 


BETR  (QEG) 


Comparisons  of  the  results  from  the  constitutive  property 
cases  PI  to  P5  leads  to  the  following  conclusions: 

1.  Increasing  (or  decreasing)  the  cohesion,  , 
increases  (or  decreases)  the  minimun  loading 
stress  for  joint  activation  while  leaving 
unchanged  che  tunnel  angles  (8)  where  joint 
activation  is  most  likely. 

2.  Reducing  the  fiction  angle  to  =  0,  thereby 
simulating  a  von  Mises  limiting  shear  strength 
of  tq/  causes  an  expansion  of  the  possible  joint 
activation  region. 

3.  Varying  Poisson's  Ratio  for  the  earth  media  from 
v  =  .25  to  v  =  .40  causes  an  expansion  of  the 
possible  joint  activation  region. 
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SECTION  3 


DYNAMIC  RESPONSE  CASES 

In  the  DIABLO  HAWK  experiment,  the  stress  wave  time  history 
involves  rise  times  to  a  peak  stress  of  0.8  kbars  in  b-10  msec. 
In  Section  4,  quasi-static  techniques  are  used  to  treat  this  rise 


Time  After  Wave  Arrival  at  O.K  kbar  Stress  Kanm 


time,  which  is  quite  lenq  with  respect  to  tunnel  response.  In 
Section  3.1,  material  properties  for  all  the  cases  are  described. 
In  Section  3.2,  six  numerical  calculations  involvin'!  much  shorter 
rise  times  of  between  LOO  and  400  Msec  are  described  and  compared 


3.1  MATERIAL  PROPERTIES 


lb A  Rock  Simulant 

In  tne  small-scale  experiments,  the  jointed  media  is  16A 
rock  simulant.  The  elastic-plastic  properties  are  modeled  with 
the  following  parameters: 

Bulk  Modulus  K  =  .131  Mbar 

Shear  Modulus  C>  =  .086  Mbar 

Poisson's  Ratio  v  =  .23 

Density  =  2.0  gm/cm3 

Yield  Strength  Y  =  .22  kbar  +  1.2P,  where  P 

is  pressure.  A  Prandtl-Reuss 
non-associated  flow  rule  is 
assumed . 

Tensile  failure  is  modeled  in  most  of  the  solutions  using  a 
tensile  failure  model  which  permits  oriented  cracks  to  develop 
(Ref.  b).  A  300  psi  =  34.3  bar  tensile  strength  was  measured 
m  a  simple  tensile  test  (Ref.  7)  on  a  16A  rock  simulant  sample. 
Most  of  our  solutions  use  tins  value;  in  one  of  the  cases,  the 
tensile  strength  is  reduced  to  6.9  bars  (100  psi)  as  a  sensitivity 
study . 

In  two  of  the  cases,  frictional  forces  along  the  sliding 
joints  are  included.  In  these  situations  the  joint  constitutive 
model  is  specified  by: 

,  =  ;  +  =  ;  +  ( tan ;  )  a 

max  o  n  o  n 

w  here 

-  maximum  allowable  shear  stress  on  joint 

mu  x 
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a  =  normal  stress  on  joint 
n  J 

r  =  the  joint  cohesion 
o 

y  =  ( tan  j )  is  the  coefficient  of  friction  and  <{> 
is  the  angle  of  friction 

The  experimental  values  (Ref.  7)  are 

i  =10  psi 
o 

f  -  30° 

The  Drucker-Prager  plastic  yield  surface  assumed  for  the  16A 
rock  simulant  material  does  not  include  shear  strength  softening 
(i.o.,  reduction  of  the  yield  strength)  due  to  plastic  shear 
deformation  or  tensile  cracking.  This  may  be  significant  in  this 
application  because  the  uniaxial  strain  loading  patii  for  this 
material  is  close  to  the  yield  surface.  Therefore,  small  changes 
in  the  failure  surface  may  lead  to  large  changes  m  the  plastic 
response  of  the  IGA  rock  simulant. 


Foam  Filler 

The  tunnel  is  filled  with 
with  the  following  properties: 

Bulk  Modulus 

Shear  Modulus 

Density 

Von  Mises  Yield 
Strength 


a  low  density  elastic-plastic  foam 

K  =  1.75  kbar 
G  =  2.5  kbar 
i>  =  .05  gm/cm' 

Y  =  500  psi  =  0.34  kbar 
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3.2  CASE  COMPARISONS  FOR  THE  0.8  KBAR  DYNAMICALLY  LOADED  TUNNELS 


Table  3  lists  the  joint  geometry,  tensile  limit,  rise  time, 
and  material  model  for  the  six  dynamic  loading  cases.  In  all 
cases  the  tunnel  is  filled  with  low  density  foam  with  properties 
described  in  Section  3.1. 

The  loading  is  a  plane  wave  of  0.8  kbar  magnitude.  The 
stress  increases  linearly  from  0  to  0.8  kbar  over  a  characteris¬ 
tic  rise  time  t  ^ ,  (100  usee  in  most  of  these  cases)  .  A  constant 

pressure  of  0.8  bar  is  maintained  after  x  .  The  far-field  con- 
c  r 

dition  of  the  pulse  (i.e.,  prior  to  interaction  with  the  tunnel) 
is  uniaxial  strain. 

3.2.1  Effect  of  Loading  Rise  Time 

In  the  first  three  cases,  the  rock  media  is  modeled  as  an 
elastic-plastic  media  with  no  joints  or  tensile  failure.  Cases 
1,  2,  and  3  differ  in  their  loading  rise  times  (t  =  100,  200, 
and  400  psec)  .  Cases  4,  5,  and  6  have  joints  which  are  perpe: 
dicular  to  the  wave  vector  as  indicated  in  Table  3;  in  these 
cases,  the  joint  properties  and  tensile  strength  of  the  media 
are  varied. 

Figure  14  shows  particle  velocity  (upper  part  of  figure) 
and  principal  stress  (lower  part)  fields  at  62  iisec  for  Case  1. 
Note  that  tile  velocity  and  stress  fields  have  a  plane  of  symmetry 
which  separates  the  two  fields.  The  scale  bars  for  the  velocity 
vectors  and  principal  stresses  are  shown  on  this  figure.  The 
following  sketch  defines  the  plotting  convention  for  the  principal 
stress : 


I'l  i  i : ■  i  I  I  i  ' 


I’l  i  1 1  r  P  vi  1  -I  i'i 

I  i  ■  1 1  ■ :  i  i '  1 1 

i  1 1.  'Win.  ini  >  niVipetH'n  t  s 
.in  pi >’’1  p iv  -  i  \<  l 

.  i  I’lip  rest.  I  i'll 


i  ■  L  i  i ’ r  in.  i  : ’.1 1 
a  i"  ■  ■  ■ 
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Table  3.  The  Dynamic  Loading  Cases 


'I 


<1 


Velocity  Field 


SO  tl/see 


I’arLicle  Velocity  .-in  <1  Principal  Stress  Fields  .it  hd  ..see  for  the 
l .  1  .is  t  i  e-i' 1 .  is  t  i  c  Med  i  a  Model  (Case  1,  0.8  kbar  l.oadiin;  in  100  ,  se 


no  loading  wave  approaches  the  tunnel  from  Hie  left.  By  02  psec, 
ho  tunnel  is  responding  very  non-symme t r ical  1  y  .  Note  that  one 
r incisal  stress  is  nearly  aero  near  the  tunnel  boundary ,  and  the 
anaont  ial  stress  is  largest  near  the  spr  i  m;ii  no  of  tire  tunnel. 

1'iuare  14  otiows  that  one  of  the  principal  stresses  is 
ii  mtly  tensile  near  Inc  tunnel  crown  (left  side  of  tunnel) . 
inure  L  i  snows  the  tensile  stress  history  at  the  crown  for 

a  son  1,  2,  and  i .  The  tensile  response  near  t:ie  crown  is  clearly 

e ;  onei  en  t  «>:i  the  stress  wave  rise  t  i  for  tiiv.es  of  a  few  hundred 
icrosecoads  .mu  less;  a  sound  wave  crosses  the  character istic 
annex  diameter  l.  n  about  22  psee .  In  fa-.  p.i.isi -static  solution 
itnout  ,  o  i ;  1 1.  activation  (see  Sect  ion  4  )  ,  t  nor-,  was  no  tensile 
tress  at  lac  crown . 

.....  ..ft  c  t  e  t  i\x  ■  c ; .  .  e  n  s  .  lx;  .  i  l  <  n  '.in  jin  .  1  o  i  n  t  shear  S  t  r  e  n  a  t  it 

fit'  WAV!.  - ..  e>  he  am  i:.x  ca-mcit.  to  all-. ■  w  tne  formation  of 
cnsiie  er.tcks,  to  c  t  .  c  ,  »t  o  tin  w  nail  • .'  f  the  cracks ,  and  to 
•a  ;  ast  stresses  an  j.icent  t  >  t  :n.  mace'  min.;  elastic  stress-strain 
e  L  a  t  i  onsn  i :  w  .  Aft.ci  a  c«.  i  I  a..;  cr  ic.wii,  :n  >  tons!  le  stress  across 
rie  crack  is  perni  t  t.od ,  a  ns  m «  smar  .'.tresses  are  permitted  on 
no  c!  .le.-.  i  •  the  ei  ick  in.  open  { tile  v.  l.itii  at  the  crack  is  con- 
inaossly  men  t  f  area)  .  ,’hosc  stress  a.i  i  uv.tnciit  s  .ire  i  r.i;  -or  tan t. 

-  'C.msc  i  c  las t  r e . pue  i i  1 1  y  .  t  1  t  c  r  s  tin  i  -  '  -c  a  1.  s  t  lens  t  i  x 1 .1  u  in  suck 
w  «y  a.,  t  )  eiiiiatice  its  own  :r>  'Wth;  thus,  my  real  i  s  t  i  c  method 

:  :  n  .act  m  :  ci.ic.-.  urowti.  must,  cons i dot  fine  a  1  t*.  red  stress  field, 

in.,  -.ami:  it  u-:i  il  .'ell  is  1 1  lowed  to  di.-vcio;  ap  to  two  cracks, 

t  i  t  . .  i  r  a  nr  <ck  t.’i'i'is,  tne  cc  1  i  is  cons idercu  to  lie  completely 
a  itt  -  i  ■  ■  .uni  no  tens  ion  at  al  i  is  tiierea  f  ter  al  1  owed .  Vne 

i  a  i  t  t'e-i  t.neiC  -  :  tens  i  le  t.iil.ire  and  tne  ;  lottiiiu  convent  ions 
.....  to  mule  t-n  the  cniCKe.i 


c<  ■  1  Is  ni  t  tie  com:  'U t a  t  i onn  1  ci a-i 


critical  tensile  1'ailure  criterion 
in  the  indicated  principal  stress 
direction ,  > 

Plane  which  experiences  tensile 
fa  1 lure . 


Crack  opens  and  relieves  the 
tensile  stress  such  that  no  forces 
or  tractions  exist  on  the  crack 
surfaces;  i.e.,  zero  normal  and 
shear  stress  on  the  crack  surface. 


In  the  directions  perpendicular 
to  the  crack  surface,  the  material 
can  maintain  tensile  stresses 
until  the  critical  tensile  failure 
criterion  i.s  exceeded  and  a  new 
crack  is  formed. 


P  1  ok  t.  i  n< !  Conven  t. ion  : 

le  crack  in  x-y  plane.  The  line  is  plotted  in 
direction  of  the  crack  (normal  to  the  principal 
ss  cans  ina  the  crack)  . 


id-plane  ctack;  i  .«• .  ,  ^  _  caused  the  crack. 


'-plane  and  x-y  plane  et  ick . 


x-y  plane  cracks ;  only  ou'-oi-pl  me  tension  is 
i  t  t  ed  . 


1  '  t  e  1  y  t  i  not  ui  ed  -  no  t  oils  l  on  pern  l  t  Led  . 


In. 


I  t  «  .it  til  .a  Ian  s  i  1 1  ■  i  imi  ie.  in  WAV!'  -  I .  (  I'Pe  . 


In  Cases  4  and  3  the 


of  the  solution  to  varia¬ 


tions  in  the  tensile  strength  of  the  rock  medium  is  examined, 
in  Case  4,  the  16A  rock  simulant  has  a  tensile  strength  of 
b.9  bars  (=  100  psi) ;  in  Case  5,  a  strength  of  34.5  bars 
1  =  500  psi)  is  used.  The  effect  on  the  tensile  crack  formation 
is  shown  in  Fiqures  17  and  18  which  show  these  two  cases  at  two 
different  times. 

Case  fa  implements  all  the  material  and  joint  properties  used 
in  tne  quasi-static  DIABLO  HAWK  simulation  (Section  4),  including 
friction  alonq  the  joints.  The  presence  of  the  friction  is  the 
only  difference  between  Case  5  and  6.  The  final  crack  pattern  in 
Case  6  is  compared  to  Case  5  in  1’iqure  19.  Fi- '  '  •.  •  -  ;  •> 

Tensile  crackiiui  of  the  material,  near  the  tunnel  surface  is 
closely  related  to  the  relative  slide  displacements  across  the 
joints.  When  slipping  occurs  across  joints,  the  plates  o!  rock 
separated  by  tne  joints  can  act  as  unsupported  beams  due  to  the 
presence  of  the  unlined  tunnel.  The  resulting  beam  bendinq  1  ends 
to  tensile  stresses  ana  fracture. 


‘\ 
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»  V  U'H  s  l  U*  M  IVU  i'.l  il 


>i  li'O  j  *  .  i  .m,l  >i)D  psi  n'.isrs  4  .in*i  n.S  kb.u  ll.b  ks  i 
■  I  is ■  * ■  •  !  «'  us  i  nr  i  n  UM  ..-.i-s  )  . 


Tensile  Strength  =  '34.5  bars  =  500  psi 


lensile  Crack  Coni  igurat  ion  at  t  =  1 2 )  esee  I  or  (ensile  Strength 

el  loo  psi  and  500  psi  (Cases  4  and  >;  O.S  khar  =  M  .  h  Its  i 
I  less  holding  in  100  ..see). 
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SECTION  4 


DIABLO  HAWK  SMALL-SCALE  STRUCTURE  PREDICTION  CASE 

This  numerical  simulation  was  intended  to  match  a  fielded 
small-scale  SRI  model.  The  calculation  results  were  obtained 
to  the  excavation  of  the  model.  The  particular  model  is 
designated  as  S-Dl-2-5  by  SRI . 

4.1  PROBLEM  CONFIGURATION 

The  tunnel  diameter  is  j" .  Joint  spacing  is  nominally 
(0.47")  .  The  boundaries  of  the  numerical  grid  are  at  9.45"  from 
tile  tunnel  center  (0 .  J  tines  the  tunnel  radius),  which  is  suf¬ 
ficient  to  minimice  edge  effects  in  the  region  of  interest.  The 
joint  geometry  is  shown  with  tiie  initial  computational  grid  in 
Figure  20. 

Material  properties  of  the  lb A  rock  simulant  and  foam  filler 
.ire  described  in  Section  3.1.  The  frictional  properties  ut  the 
joints  are  also  discussed  m  that  section. 

The  tunnel  and  jointed  media  are  subjected  to  a  planar  stress 
wave  of  0.8  kbar  peak  amplitude.  The  rise  time  is  5-i.O  ns.-c 
and  decay  10  0  :;;s.,c .  Tnese  times  are  very  lone  compared  to 
tiie  22  ..sec  sound  wave  transit  time  across  tile  tunnel  diameter 
in  tiie  IGA  material.  Tiie  actual  tunnel,  therefore,  is  in  a 
"guas l- static"  state. 

4.2  QUASI -STATIC  SOLUTION  TECHNLQFF 

Tiie  juasi-static  solution  is  found  as  indicated  on  Figure  21: 

1.  Obtain  the  eLustic  solution  to  a  tunnel  loaded  witii 
tiie  0.8  kbar  (far-tield  uniaxial  strain)  stresses. 


Figure  21.  Successive  Steps  in  Quasi-Static  Numerical  Simulation  Technique. 


2.  Using  tne  elastic  solution  as  initial  conditions, 
activate  the  elastic-plastic  material  model  and 
calculate  the  dynamics  until  a  steady  state 
solution  is  obtained. 

3.  Using  this  elastic-plastic  solution  as  initial 
conditions,  a  calculation  is  performed  with  the 
slide  lines  introduced  to  model  the  joints  and 
permit  block  motion  and  deformation.  Also,  the 
material  may  fail  in  tension  if  the  tensile 
stresses  exceed  500  psi  in  any  principal 
direction.  The  dynamic  calculation  is  continued 
until  tensile  failure,  block  motion  and  tunnel 
deformation  have  stopped,  and  a  steady-state 
solution  is  thereby  obtained. 

4.3  JOINT  SLIDING 

Peak  slide  displacements  of  l'i  of  the  tunnel  diameter  are 
predicted  in  the  quasi-static  solution.  In  a  full  sized  tunnel, 
with  a  diameter  of  10  feet,  this  corresponds  to  1.2  inches  of 
sliding  . 

Figure  22  shows  the  slide  displacement  along  the  joint  as  a 
function  of  distance  (scaled  by  tunnel  radius)  from  the  tunnel 
axis.  The  sliding  is  a  maximum  at  about  one  radius  along  the 
joints  nearest  the  tunnel  crown  and  invert  (Joints  3  and  9  in 
Figure  22) . 

4.4  TENS ILL  CRACKING 

The  tensile  cracking  is  closely  related  to  the  relative  slide 
displacements  across  the  joints.  When  slipping  occurs  across  a 
joint,  the  plates  of  rock  separated  by  the  joints  act  as  beams. 
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The  presence  of  the  unlined  tunnel  removes  localized  support, 
and  beam  bending  leads  to  tensile  stresses  and  fracture. 

The  crack  pattern  obtained  in  the  numerical  solution  is 
shown  in  Figure  23.  The  "beams"  terminating  in  the  opening 
show  significant  cracking  of  the  unsupported  cantilevered  ends. 

4.5  TUNNEL  CLOSURE  AND  COMPARISON  WITH  EXPERIMENTAL  DATA 

The  tunnel  closure  predictions  arc  shown  on  Figure  24.  The 
closure  along  the  vertical  axis  (i.e.,  crown-invert)  is  approxi¬ 
mately  0.55.  The  maximum  closure  of  1.2%  is  located  at  0  =  30°. 

Figure  24  shows  a  comparison  of  the  actual  and  predicted 
tunnel  closure  (as  a  %  of  the  diameter)  versus  tunnel  angle.  The 
agreement  is  quite  good  except  near  the  crown-invert  (0  =  0°  and 
180°).  The  experimental  data  show  a  1.6%  closure  at  the  crown- 
invert  compared  to  a  0.5%  predicted  value.  This  difference  is 
probably  due  to  the  moderate  sliding  on  the  joints  and  associated 
tensile  failure  and  relaxation  which  occurs  in  the  calculations, 
but  not  in  the  experiment. 

In  the  quasi-static  calculation,  tensile  cracking  occurs 
after  the  joints  are  "unlocked"  and  joint  sliding  begins  (see 
Section  4.4  and  Figure  22).  Tensile  failure  is  not  evident  in 
the  post-shot  experimental  model  (Figure  25).  If  the  metal  casing 
(which  is  not  modeled  in  the  calculation)  reduced  the  sliding  on 
the  joints,  or  if  the  actual  joints  are  stronger  (higher  cohesion 
and/or  higher  coefficient  of  friction)  than  used  in  the  calcula¬ 
tion,  then  there  would  be  less  block  motion  in  the  experiment, 
and  therefore  there  would  also  be  less  tensile  failure. 

Figure  25  shows  the  post-shot  experimental  model.  Note  the 
deformation  of  the  thin  (.003  inch)  metal  casing  (or  liner) 
surround  i.  nu  the  foam.  This  casing  may  have  inhibited  sliding 
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small-scale  model  ot  tile  amount  of  relative  displacement  across 
tiie  joints  is  iiiqhly  desirable  for  comparisons  with  the  WAVh-L 
code  predictions  of  Fiqure  2.1. 
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APPENDIX  A 

JOINT  ACTIVATION  CONDITIONS  NEAR  A  CYLINDRICAL  TUNNEL 


The  relationships  between  the  principal  stress  components 
(ni,  o 2 ) ,  joint  and  tunnel  angles  (  t,  B ) ,  and  constitutive  acti¬ 
vation  properties  will  be  obtained  using  a  Mohr  diagram  (Ref.  8). 
The  joint  constitutive  model  is  specified  by 


max 


i  + 

o 


no 


n 


or 


where 


max 


+  (tan  ■+• )  o 


(A- 1 ) 


=  the  maximum  allowable  shear  stress  on  the  joint 
x  J 

=  the  normal  stress  on  the  joint 
n  J 

=  the  joint  cohesion 
o 

n  =  tan  ;  is  the  coefficient  of  friction  and  :  is 
the  angle  of  friction. 


Figure  A- 1  (page  57)  is  a  Mohr  diagram  which  shows  (1)  the 
shear  stress  and  normal  stress  ■'  for  the  joint  constitutive 
model  (Equation  A-l),  and  (2)  the  possible  shear  and  normal  stress 
states  corresponding  to  the  principal  stresses  i  and  o .  .  The 
following  sketch  (page  54)  indicates  the  geometrical  interpretation 
of  the  Mohr  diagram  in  terms  of  the  stresses  on  a  joint  near  a 
tunnel  surface: 
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Doing  tho  algebra  and  solving  for  a? ,  we  find, 


1 


/ sin(2y-;)  +  sin  A 
;  \ sin  (2 y-f)  -  sm  v 


o  / _ sin  _ 

n  (sin  ( 2  y  - ;  )  -  sin  <;. 


Equation  A-3  relates  the  two  principal  stresses  such  that 
an  oriented  joint  (with  angle  y  between  the  tunnel  normal  and 
the  joint)  will  lead  to  joint  activation.  A  more  physical  and 
useful  relationship  for  a  plane  stress  wave  engulfing  a  tunnel 
can  be  obtained  on  a  <>:»:•  '  a  a  • 

On  an  .  tunnel  surface,  one  principal  stress  component 

vanishes;  i.e.  ,  Ji  =0.  Thus,  the  tangential  stress,  u2  r  on  the 
tunnel  surface  for  joint  activation  is 


\  2 1  cos  ; 

_ _  _ sin:  )  _  _ o _ | _ 

s  i  n  ( 2  ,  -  ;  )  -  sin  ;  sin(2y~v)  -  sin 


(A-4) 


When  a  plane  wave  engulfs  a  tunnel,  however,  the  tangential  sc 
varies  with  the  angle  on  the  tunnel  surface  as  illustrated  in  the 
following  sketch: 


=  A ( .■ )  ' 


[ 1  +  2cos (2  ■  )  | 
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•nui  Lo.niing  Stress.  '  on  Tunnel  Surlaee 
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In  the  sketch,  A(8)  is  the  amplification  factor  of  the 
°load  wave •  A(8)  depends  on  the  angle  B  on  the  tunnel  surface 

and  on  the  elastic  Poisson's  Ratio  (v)  of  the  media.  The  deri¬ 
vation  for  A(B)  can  be  obtained  by  superposition  of  elastic 
solutions  derived  in  Reference  8.  Figure  5  (in  body  of  report) 
shows  A (8)  for  v  =  .25  and  .4. 

Since  o  i  =0  and  o2  =  A(8)  ,  °^oacj'  Equation  A-4  shows  that 


load 


2t  cos  di 

_  U  2  _  _ _ o _ _ _ 

a(B)  [sin  (2tp-<p )  -  sin  $]A(f3) 


where 


A(B)  =  2 


.  1  +  2cos  (28)) 


l-2v 


This  is  the  activation  stress  load  relationship  (Equation  1) 
used  in  the  report. 


(compress ive ) 
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